EXISTENCE OF TRAVELLING WAVES IN DISCRETE 
SINE-GORDON RINGS 



GUY KATRIEL 



Abstract. We prove existence results for travelling waves in discrete, 
damped, dc-driven sine-Gordon equations with periodic boundary con- 
ditions. Methods of nonlinear functional analysis are employed. Some 
unresolved questions are raised. 



1. Introduction 



The damped, dc-driven discrete sine-Gordon equation, known also as the 
driven Frenkel-Kontorova model, with periodic boundary conditions, arises 
as a model of many physical systems, including circular arrays of Josephson 
junctions, the motions of dislocations in a crystal, the adsorbate layer on the 
surface of a crystal, ionic conductors, glassy materials, charge-density wave 
transport, sliding friction, as well as the mechanical interpretation as a model 
for a ring of pendula coupled by torsional springs (we refer to [SJ ITT] 
and references therein). This model has thus become a fundamental one for 
nonlinear physics, and has been the subject of many theoretical, numerical 
and experimental studies. The system of equations is 

(1) tjJt + Tfy+wa&j) =F + K[4> j+l -24> j +<t> i - l }, V? e Z 

with the parameters r > 0, K > 0, F > 0, with the periodic boundary- 
condition 



where m > 1 (we note that in view of the boundary conditions we are re- 
ally dealing with an n-dimensional system of ODE's rather than an infinite- 
dimensional one). In numerical simulations, as well as in experimental work 
on systems modelled by (Q,©, it is observed that solutions often converge 
to a travelling wave: a solution satisfying 
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where the waveform / : K — ► R is a function satisfying 

(4) f(t + T) = f(t) + 2tt VteK. 
The velocity of the travelling wave is given by 

(5) v=jr- 

However, as has been pointed out in |10| . even the existence of such a solution, 
has not been proven, except for the case of small K in which existence of a 
travelling wave for some values of F had been proven in 0] . 

In the 'super-damped' case, in which the second-derivative term in Q is 
removed, there are very satisfactory results about existence and also global 
stability of travelling waves (p^, theorem 2). Such results rely strongly on 
monotonicity arguments. Recently Baesens and Mackay [2] have managed to 
extend these arguments to the 'overdamped' case of their results apply 
when 



(6) r > 2V2ir+T, 

and say that there exists a travelling-wave solution which is globally stable if 
and only if Q,© does not have stationary solutions. We do not know whether 
in general the non-existence of stationary solutions implies the existence of a 
travelling wave. 

We note that a function / is a waveform if and only if it satisfies 10} and 

(7) /"(*) + Tf'(t) + sin(/(t)) =F + K[f(t+^T)- 2f(t) +f(t- ^r) 

Here we obtain several existence results for travelling-waves under condi- 
tions not covered by the existing work, described above. 

Theorem 1. Fixing any T > 0, K > 0, and given any velocity v > 0, there 
exists a travelling-wave solution of PJ), @) with velocity v for an appropriate 
F > 0. 

Theorem 2. For any F > 1 there exists a travelling-wave solution of f7)) . 

m. 

Theorem 3. Assume that n does not divide m. Fixing any F > and T > 0, 
for all K sufficiently large there exists a travelling-wave solution of 0), 0) 
for any F>F,T> f. 

Theorem 4. Fixing any F > and K > 0, for all V > sufficiently small 
there exists a travelling-wave solution of 0), 0) for any F > F , < K < K . 

We remark that the assumption that n does not divide m cannot be re- 
moved from theorem |3 since if n divides m the coupling term vanishes and 
0,0) reduce to the equation of a running solution of a dc-forced pendulum, 
which, fixing T > 0, is known to have a solution only when F exceeds a 
positive critical value [S]- 
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It is interesting to note that theorem 01 demonstrates that for some param- 
eter ranges there is coexistence of stationary solutions and travelling waves of 
0, ©■ Indeed, it is well known [3] that, fixing K, for F sufficiently small 
there exist stationary solutions of (JTJ , 10) , and these obviously do not depend 
on r. Hence we can take T > sufficiently small so that theorem 0] ensures 
also the existence of travelling waves. This phenomenon cannot happen in 
the super-damped case, nor in the overdamped case in which © holds, since 
in these cases existence of stationary solutions implies that the w-limit set of 
every orbit is contained in the set of stationary solutions ^ |2] • 

Along the way we will prove that 

Proposition 5. An upper bound for the velocity v of any travelling wave is 
given by 

(8) «<-. 

and a lower bound, in the case F > 1, is given by 

F- 1 



(9) v> r . 

In the next section we prove the results stated above. In section|3we discuss 
the meaning of our results in connection with existing numerical studies of 
the discrete sine-Gordon equation, and point out some further mathematical 
questions which arise from our results and remain open. 

2. PROOFS OF THE RESULTS 

Our method of proof involves re-formulating the problem as a fixed point 
problem in a Banach space, and applying results of nonlinear functional anal- 
ysis. Our approach is thus close in spirit to 6 , which deals with travelling 
waves in globally coupled Josephson junctions. 

We transform the problem (0J,I7J) by setting 

f(t) = u(vt)+vt, 
where the wave- velocity v is defined by © and u satisfies: 

(10) u(z + 2tt) = u(z) VzeR. 
JZJ) can then be written as 

v 2 u"(z) + Tvu'{z) + sin(z + u(z)) 

(11) = F-Tv + K u(z + 2ir— ) -2u(z) + u(z-2tt— ) 

L V n J \ n ) 

Dividing by v 2 and setting 





„ m\ 


u^z 


+ 2tt— 




n ) 




1 


A 






V 
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we re- write 1(11(1 in the form 

u"(z) + XTu'(z) + X 2 sm(z + u(z)) 

r / Tfl \ f TTX 

(12) = \ 2 F- \T + \ 2 K u(z + 2it— )-2u(z) + u \z-2-k— 

-V n / V n 

We note that if u(z) satisfies (|10ll . (ll2(l then so does u(z) = u(z + c) + c, for 
any c£l. Thus by adjusting c we may assume that u satisfies 

(13) / u(s)ds = 0. 



We note now that if u satisfies <|10|) . <|12ll . then by integrating both sides of 
(|12ll over [0, 2ir] we obtain 

(14) F=^ + ^- [ sm(s + u(s))ds. 

A 2n J Q 

We can thus re-write lfi"2"|l as 

1 r 27r 

(15) u"(z) + XTu'iz) + A 2 sin(z + u(z)) = A 2 — / sin(s + u{s))ds 

+ \ 2 K uiz + 2-K— )-2u(z) + u [z-2ir— 
-V n ) \ n 

Conversely, if u satisfies (|14H and l|15|) then it satisfies (|12|) . We have thus 

reformulated our problem as: find solutions (A,m) of H1U|) . I|13|I . (|14|I . (|15|) . The 

idea now is to consider A as a parameter in (|15|l and try to find solutions u 

satisfying (|10|I . 113|I . 115(1 and then substitute A and u into 1)14(1 to obtain the 

corresponding value of F . This is the same idea as used in the numerical 

method presented in [jy, but here it is used as part of existence proofs. We 

claim that 

Proposition 6. For any value \, there exists a solution u of \T$j) satisfying 

We note that this proposition immediately implies theorem ^ since given 
any v > it shows that we can solve ((15|) with A = ^ , hence obtain a travelling 
wave with velocity v, for the value of F given by ((1411 . 

To prove proposition [S] we will use the Schauder fixed-point theorem. We 
denote by X, Y the Banach spaces of real- valued functions 

X = {ue H 2 [0, 2tt] I u(0) = u(2tt), u'(Q) = u'(2ir), / u(s)ds = 0}, 

u(s)ds = 0}, 



with the norm 



Y = {ue L 2 [0,2tt] I f 
Jo 

\\u\\ Y =(^J 2 J {u{s)) 2 ds)\ 
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and by L\ : X — ► Y the linear mapping 

L\(u)(z) = u"{z) + XTu'(z) - X 2 K \u[z + 2tt — J - 2u(z) + u[z - 2tt— J 

We want to show that this mapping is invertible and derive an upper bound 
for the norm of its inverse. Noting that any u S X can be decomposed in a 
Fourier series u(z) = a i e%lz (with a_; = af), we apply La to the Fourier 

elements, obtaining, 

L X\ e ) — Hl e > 

where 

W = -I 2 - 2KX 2 ( cos f — ) - l) + MTi, 



so that 

(16) | w |= [(> + 2XA 2 (cos(^) -l)) 2 + AW 

which does not vanish if T > 0. Thus the mapping L\ has an inverse satisfying 
L^ 1 (e llz ) = -j^e Uz ■ Since L^ 1 takes y onto X, and since X is compactly 
embedded in Y, we may consider L^ 1 as a mapping from y to itself, in which 
case it is a compact mapping. We also note using (O that 

(17) ^W< ? ax^<^. 

We also define the nonlinear operator N : Y — > y by 

1 f 2 " 

N(u)(z) = — sin(z + it(z)) H / sin(s + u(s))ds. 

27r Jo 

It is easy to see that N is continuous, and that the range of N is contained 
in a bounded ball in Y, indeed we have 

/ 1 f 27T \ 5 

|| sin(z + u(z))|| Za = {^J ( sin ( s + w(s))) 2 ds) 2 < 1, 

and since N(u) is the orthogonal projection of — sin(z + u(z)) into Y, we have 

(18) \\n{u)\\ y < i Vu e y. 

We can now rewrite the problem H1U|) . I|13[) . (|15I) as the fixed-point problem: 

(19) u = X 2 L^ 1 oN{u). 

The operator on the right-hand side is compact by the compactness of , 
and has a bounded range by H17|) . l|18fl . so t na ^ Schauder's fixed-point theorem 
implies that l|19[) has a solution, proving proposition [5] (we note that by a 
simple bootstrap argument a solution in Y is in fact smooth). Moreover, 
defining 

£ = {(A, u) € [0, oo) x y | u = X 2 L^ o AT( U )}, 
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Rabinowitz's continuation theorem |Jj implies that the connected component 
of X containing (A, u) = (0, 0), which we denote by C, is unbounded in [0, oo) x 
Y. Since for any A > we have, from (JTSJl, jT^I , the bound < ^ for 

solutions (A, u) of l|19|l with A S [0, Ao], the unboundedness of the set C must 
be in the A-direction, that is, there exist (A, u) € C with arbitrarily large 
values of A. 

We can now consider the right-hand side of (fTIjl as a functional on [0, oo) x 

Y: 

r i f 27r 

(20) $(A,u) = - + — / sin(s + u(s))ds, 

A 27r ,/ 

and our strategy in proving theorems 12141 is to prove solvability of the equation 



(21) $(A,tt)=F, (A,ti)6S 

(in fact we shall prove solvability of H21|) with X replaced by C C X). We note 
that by the boundedness of the sine function we have 

(22) lim $(A,it) = +oo, 

A^0+, (A,u)6C 



(23) limsup $(A,it)<l. 

A^+oo, (X,u)eC 

Since C is a connected set and $ is continuous, (|22l) implies that 
Proposition 7. For any F satisfying 

(24) F>F= inf $(A,w), 

— (A,u)6C 

£/iere exists a travelling wave. 

Since (|23|l implies that F < 1, this proves theorem |2j 

We now prove the lower and upper bounds for the velocities of travelling 
waves given in proposition [5] These follow from ^.1)1 4(1 and 

Proposition 8. For any (A, u) £ X with A > we have 

< - < $(A,it) < T + 1. 
A A 

The upper bound follows immediately from the definition H2()|l of <&(A, u) 
since j- J Q * sin(s + u(s))ds < 1. The lower bound follows from the claim that 



(25) 



(A, u) G X =*> / sin(.s + u(s))ds > 0. 

./0 
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To prove this claim we multiply fT5|l by 1 + u'{z) and integrate over [0, 2tt] 
noting that 

U\ S 



2tt— )ds 



2w—ju'(s)ds = J u(s)u'[ 

i (s)uys — 2tt — Jds 



so that we obtain 

sin(s + u(s))ds. 



(A,u)e^rl| ( u '( s) )2 ds = A i_^ 



This proves (|25|) since the left-hand side is non-negative and cannot vanish 
unless u = 0, but (A, 0) £ £ for A > 0. 

We now turn to the proof of theorem [21 

Proposition 9. Assume n does not divide m and V > 0. Given any Ao > 
and e > 0, there exists Kq such that for K > K , T > T we have that 

1 f 2 " 

(26) — / sin(s + u(s))ds <e i/(A ,«)eE. 
2tt Jo 

To see that proposition [5] implies theorem |3 we fix some F > 0, T > 0, 
and assume r > T. We choose An > £ and set e = F — -r-. We then choose 

— F Ao 

ifo according to proposition El so that (j26(l holds, which implies that when 
K > Ko we have $(Ao, u) < F for any u with (Ao, u) G C. Thus F_< F, where 
F_ is defined by (|24|) . so proposition[71 implies the existence of a travelling wave 
for any F > F. 

We now prove proposition El Let Ao > and e > be given. Assume 
(Ao,u) £ S, so that (|19|l holds with A = Ao. Let (m,n) denote the greatest 
common divisor of m, n and let 

m n 

P= 7 r, 9= 7 r- 

(m, nj (m,n) 

Since we assume n does not divide m we have q > 2. Let Fo be the subspace of 
y consisting of ^-periodic functions, and let Y\ be its orthogonal complement 
in Y . We denote by P the orthogonal projection of Y to Yq. Setting 

u = P{u), Ul = (I-P)(u), 

we have u = Uq + U\ with Uq € Yq, u± € Y\. Applying P and / — P to (|19|) . 
and noting that La commutes with P, we have 

(27) « = AoL A >FoiV( tlo + tll ) 

(28) Ul = X 2 L^ o (I - P) o N(u + ui) 
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We will now use l(TH)) to derive a bound for \\L Xo |y x 11^,^ which goes to as 
K — ► oo. We note that 

(29) ||^VJk,y i < & m« l] ^ 1 

so we need to find lower bounds for the |/t/|'s for which q does not divide I. 
We define 

/ 2npl \ 

p = max cos 

J>1, ?/Z \ q J 

and note that since p, q are coprime we have p < 1. 
We define 

a = 2K\l{\ - p) - VK, 

and we shall henceforth assume that K is sufficiently large so that a > 0. For 
each I > 1 we have either I 2 < a or I 2 > a, and we treat each of these cases 
separately. 

(1) In case I 2 < a, we have 

l 2 + 2KX 2 (p- 1) < -s/K, 
and by the definition of p 



/ 27rmZ \ / 2npl N 

cos = cos f 



\ n I \ q 

so that 

/ 2 + 2ifAg(cos(^-^) -l) < -VK, 
which by (|16|) implies 

(30) | > 
(2) In case / 2 > a, we have, since implies |/t;| > XoTl: 

(31) |w| >A r^> Aofv^Aof [2ATAg(l-p)- 

From (|30(I . H31| ) we obtain that hin^^oo |/t;| = +oo uniformly with respect to 
I > 1 which are not multiples of q, hence by l|29|) 

lim HL^IyJI^y^O. 

K — >oc u 

In particular we may choose K such that for K > K we will have 



ll L A Vilk,Yi < T2- 
^0 

By (ESJl and JTSJl this implies 
(32) ||«i||y<c 
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Thus 



— / sin(s + u(s))ds < — / sin(s + uo(s))ds 
27T J 2tt J 



2tt 



— / [sin(s + u(s)) — sm(s + uo(s))]ds 

Z7T 



(33) 



< 



1 

2^ 



sin(s + uo(s))ds 



1 

2^ 



\ui(s)\ds 



From (|32[) and the Cauchy-Schwartz inequality we have 



(34) 



271 



i- / ,/,(.s)|,/.S-_ 



/27T 



2 77 



(ui(s)) 2 ds <e 



From trigonometry we have 



sin(s + uo(s))ds — / sin(s) cos(uo(s))ds + / cos(s) sin(uo(s))ds, 



but the functions cos(uo(s)), sin(uo(s)) are ^-periodic with q > 2, which 
implies that they are orthogonal to cos(s), sin(s), so that we have 

f277 

sin(s + u Q (s))ds = 0, 

which together with (|33|l and 1134(1 implies (|26J) . concluding the proof of propo- 
sition El 

We now turn to the proof of theorem 0] We first note that from we 
have 

M>|z 2 + 2XA 2 (cos(^)-l), 
so that if we assume 



< A < A 



1 1 

< 



then we have \m\ > h for all I > 1, hence ||L A 1 ||y : y < 2, independently of T. 
From l|19f) we thus have 

(X,u) e E, < A < A => ||u||y < 2A 2 . 

We now choose Ai < Aq so that 2A 2 < ^F. Thus (\i,u) £ E implies that 



(35) 



1 

2* Jo 



sin(s + u(s))ds 





-f 




2* Jo 



[sin(s + u(s)) — sin(s)]ds 
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Finally, we choose To so that 

(36) ^<\f. 

(|55j) . (|Sni> thus imply that when < T < T 

(Ai,u) eS^ $(Ai,u) < F, 

so that we have F < F, where F_ is defined by (J2H, hence proposition [7| 
implies the existence of a. travelling wave for any F > F. 

3. Discussion and further questions 

In the numerical and experimental explorations of the dynamics of sine- 
Gordon rings a useful method of representation consists in dis- 
playing the velocity-force characteristic. In the case of the travelling waves 
studied here, since the velocity of the waves is given by v = 4, the velocity- 
force characteristic is the subset of the (F, u)-plane given by 



{(•(?•)■•)) I (;-M 



where the set £ and the functional <E> are as defined in the previous section. 
Examining the velocity-force characteristic as numerically computed in |Hj 
(fig. 2), we see that F is a non-monotone function of v. This means that 
for some values of F equation l|21|l has more than one solution, or in other 
words that there exist multiple travelling waves with different velocities for 
the same value of F. On the other hand, the fact that according to the 
available numerical evidence F is a function of v leads us to conjecture that, 
fixing r > 0, K > 0, for each given velocity v > there is a unique travelling 
wave with velocity v, for an appropriate F. Thus, our conjecture is that 
uniqueness holds in theorem^ or in other words that the fixed point problem 
(|19(1 always has a unique solution. Let us note that for < A < T it is easy to 
show, using Ijl7|l . that the right-hand side of 1)19(1 is a contraction from Y to 
itself, hence we may replace the use of the Schauder fixed-point theorem by 
the Banach contraction-mapping principle, which implies uniqueness. Thus, 
at least for velocities v > jt, we have uniqueness in theorem^ However for 
lower velocities this argument does not work so a proof of the above conjecture 
will require some new idea. 

The non-uniqueness of travelling waves for some values of the parameters 
r, K, F, mentioned above, implies important consequences for the dynamics 
of the discrete sine-Gordon ring, such as instability of some of the travelling 
waves, and bistability of travelling waves leading to hysteresis as the force 
F is varied. It would be interesting to determine whether these phenomena 
can occur in the large K and the small T regimes for which existence of 
a travelling wave has been proved in theorems 13141 Moreover, stationary 
solutions and travelling waves do not exhaust the dynamical repertoire of the 
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discrete sine-Gordon equations in the under-damped case: quasi-periodic and 
chaotic behavior is reported [HI El An interesting question is to determine 
conditions on the parameters which ensure that there exists at least one locally 
asymptotically stable travelling wave. 

Returning to the issue of existence of travelling waves, which has been the 
focus of our investigation, we note an intriguing question which arises from 
our results, and remains unanswered. 

Let us define, for fixed r > 0, K > 

F (r, K) = M{F > | a travelling wave of Q}, © exists}. 

Theorem implies that F (r, K) < 1 for all V > 0, K > 0. Theorem implies 
that, when n does not divide m, Vocik^oc Fq(T, K) = 0. Theorem 0] implies 
that limr-»o Fq(T, K) = 0. Is it true, though, that for each r > 0, K > we 
have Fo(T, K) > 0? In other words, is it always true that (fixing r > 0, K > 0) 
for sufficiently small F > a travelling wave does not exist? We have not 
been able to prove or disprove this conjecture, and can only offer the following 
remarks: 

(i) If T, K satisfy ©, then indeed -F (r, K) > 0, since for sufficiently small F 
there exists a stationary solution of Q),©, so the results of |2J imply that no 
travelling wave exists for small F. However, as we have remarked, theorem 0] 
shows that in general travelling waves and stationary solutions may coexist. 

(ii) If n divides m then, as was noted in the introduction, the existence of 
travelling waves reduces to that of running solutions of the forced pendulum, 
hence it is well known that Fq(T,K) > for all T,K. However, in the case 
that n divides m we also have noted also that the result of theorem |31 does 
not hold, so that the case that n divides m is rather special and may not be 
indicative of the general case. 

(iii) The conjecture that Fq(T, K) > is supported by the notion of 'pinning' - 
the phenomenon whereby travelling waves are unable to propagate in discrete 
systems when the applied force is small. However, whether this effect indeed 
holds in general in underdamped systems (as opposed to the overdamped 
case - see (i) above) is unclear to the best of our knowledge. Moreover, it 
is conceivable that Fo(T, K) = but pinning still occurs - if for small F a 
travelling wave exists but is unstable. 

(iv) Since, by proposition |H1 we have $(A, u) > for all (A, u) e E, A > 0, we 
have that Fq(T, K) = if and only if 

liminf <&{\,u) = 0. 

A^ + oo, (A,m)GS 

Thus, determining whether the above equality can hold could be a route to 
resolving our question, but we have not been able to do so. 
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We conclude with one more question: clarify the connection, if any, between 
the travelling waves obtained in 0] for small values of K > and those 
obtained by us for large values of K in theorem [3J 
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